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■ Abstract 

(N 

Using a regularization by putting the system in finite volume, we develop a novel 
approach to form factor perturbation theory for non-integrable models described as per- 
turbations of integrable ones. This permits to go beyond first order in form factor pertur- 
bation theory and in principle works to any order. The procedure is carried out in detail 
for double sine-Gordon theory, where the vacuum energy density and breather mass cor- 
rection is evaluated at second order. The results agree with those obtained from the 
truncated conformal space approach. The regularization procedure can also be used to 
compute other spectral sums involving disconnected pieces of form factors such as those 
• that occur e.g. in finite temperature correlators. 

^ ■ 1 Introduction 

> 

Q\ | Form factor perturbation theory (FFPT) was developed in [T] in order to evaluate quanti- 

ties in a non-integrable model obtained as a perturbation of an integrable one. Writing the 
, Hamiltonian in the form 

o 

ON 

o 



-^nonintegrable ^integrable + ^ f dx^f(t, : 



where ^ denotes the local perturbing field that breaks integr ability, the first order corrections 
to the vacuum (bulk) energy density and particle masses are given as 

V- 

5£ vac = A(0|*|0) A=0 
5M 2 ab = 2AF^(ivr,0) 

where 

Ft. in • • • M = <0|*(0, 0)\A h (<?i) ...A ln (tf n )) A=0 (1.1) 

are the form factors of the perturbing operator calculated at the integrable point A = and b 
denotes the charge conjugate of particle species b. It is possible to evaluate first order correc- 
tions to the two-particle S matrix and also the widths of decays induced by the perturbation 

The evaluation of higher order corrections has not been developed; simple considerations 
along the lines of [Tj lead to divergent expressions. However there is no place for mass renor- 
malization by counter terms analogous to standard Feynman perturbation theory because the 
operator (jl.ip defined by the form factors is already well-defined and physical. This is also 
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confirmed by the fact that when the non-integrable model is formulated using the truncated 
conformal space approach pioneered by Yurov and Zamolodchikov [3] the mass gaps turn 
out to be finite and well-defined (the vacuum energy can still have divergent contributions 
depending on the ultraviolet weight of but the differences between energy levels are all 
finite). 

This leads to the central idea of the paper: since the TCSA expression for the relative 
energy levels in a finite volume L is finite, and the ingredients necessary to evaluate finite 
volume perturbation theory can be determined from infinite volume form factors using the 
approach developed in [HE], one can write down finite and well-defined analytic expressions 
for the perturbation of finite volume energy levels (which are accurate up to so-called residual 
finite size effects that decay exponentially with the volume i.e. are non-analytic in 1/L, 
i.e. valid to all orders in 1/L). Then the quantity of interest (bulk energy density or mass 
correction) can be expressed directly in finite volume and the infinite volume limit is taken 
only at the end of the calculation. This is the same philosophy that was used to obtain the 
expression of finite temperature one-point functions in [5]. 

Eventually, since first order FFPT was used in [5] to derive the expressions of one-particle 
and two-particle diagonal matrix elements in finite volume, nothing new is to be gained from 
the application of finite volume techniques at first order. However, we get new results at 
second order: a consistent, generally valid way of calculating corrections to vacuum energy 
density and particle masses. It can also be extended to other quantities such as the S matrix, 
and to higher order FFPT corrections as well. 

It is best to consider a concrete model to develop and test this approach. The model of 
choice is the double sine-Gordon model defined by the Hamiltonian 

Hbsg = j dx Q (d t (p) 2 + i (d x (p) 2 - [i : cos (3ip : +A : sin :^ (1.2) 

understood as a perturbation of the massless free boson (which also defines the normal or- 
dering). It has attracted interest recently chiefly because it is a prototype of non-integrable 
field theory which can be understood by application of techniques developed in the context of 
integrable field theories [6j [9] and it also has several interesting applications [U EJ [8] such as 
to the study of massive Schwinger model (two-dimensional quantum electrodynamics) and a 
generalized Ashkin- Teller model (a quantum spin system) which are discussed in [6]. 

The double sine-Gordon model (II. 2|) can be considered as a non-integrable perturbation of 
the integrable sine-Gordon field theory obtained by setting A = [6j. Form factor perturbation 
theory was applied to the double sine-Gordon model in [9J; for the particular version in eqn. 
(|1.2|) it was shown that the corrections to the breather masses vanish to first order in A. 
However later semiclassical considerations [10] seemed to contradict these naive expectations, 
yielding mass corrections which were of first order in the coupling A. In [12] it was shown that 
a precise numerical determination of the spectrum contradicts this conclusion and upholds the 
naive picture obtained from form factor perturbation theory: i.e. there are only second order 
corrections, and in fact all odd orders vanish since the entire spectrum turns out to be even 
under A — ► —A. However, at that time the mass correction could not be calculated theoretically 
due to the lack of FFPT beyond first order. Therefore this model is an interesting testing 
ground for the present work, and it is also made ideal by the absence of first order corrections 
which makes comparison to numerical results easier. The numerical results which are compared 
with the theoretical predictions are obtained from TCSA which was first developed for the 
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sine-Gordon model in [TT] and generalized to the double sine-Gordon model in [9]; to achieve a 
better precision we use an improved version developed for the work [12J and described therein. 



2 Bulk energy correction 

The general formula for second order corrections to energy levels can be written as 



<*i#iifc>r 



E. 



(o) 



E 



(0) 



where E-f^ are the unperturbed energy eigenvalue corresponding to the eigenstate \i) and H\ 
is the perturbation to the Hamiltonian. In our case 



f L 
Hi = A / dx : sin —ip : 

Jo 2 

where L is the spatial volume of the system. With periodic boundary conditions the matrix 
elements of Hi vanish unless the momenta of states \i) and \k) coincide; therefore when i is 
taken to be the vacuum, only states with zero total momentum contribute. In addition, the 
topological charge of \k) must also be zero, otherwise the amplitude vanishes. Furthermore, 
Hi is odd under C : ip — > —ip and so the C-parity of the contributing state must be odd 
as well (the nth breather B n has C-parity (—1)™). Only contributions of breather states are 
necessary to evaluate because our numerical data will come from a part of the attractive 
regime of sine-Gordon theory where solitons are heavy and contribute little to the summation 
over k, well below the available numerical precision. 
When the momentum of the state \k) is zero, 



(0|JZi|fc) = A(0 
is independent of x, therefore 



k 



SE n 







: expi^ip(0, 0) : 



#(°) _ 



(2.1] 



where the subscript L designates finite volume matrix elements (the two exponential terms 
in the sine give equal contributions due to parity). For £ < 1/20, the lowest lying states 
contributing to the sum are 

1^(0)), |B 3 (0)), 

\Bi(9i)B 2 {02)) with mi sinh B\ + m 2 sinh 9 2 = 
and \Bi{6i)Bi(p<2)Bi{Qz)) with mi sinh 6\ + mi sinh 6 2 + mi sinh #3 



where 



^or 







TOfc = 2Msin 



the notations M and ^ cf. Appendix [XI 
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are the breather masses, and the rapidities of the breathers are indicated in parentheses. In 
details 



SE (L) 



-X 2 L 2 



Oi 

0i fii 





: expi^(/?(0, 0) : 


ft(0)) L 


2 




: expi^y?(0, 0) : 


bm) l 








\ 2 L 2 - 









mi 



m 3 



:expif^(0,0) : B 1 {6 1 )B 2 {e 2 ) 



(mi cosh 9\ + ni2 cosh 02 ) 
:expifvp(0,0) : B 1 (9 1 )B 1 (9 2 )B 1 (9 3 ) 



(nil cosh 6*i + nil cosh 02 + f»i cosh #3) 



+ O (e"^) + . . . 



where the presence of correction terms decaying exponentially with the volume is indicated, 
and the ellipsis denotes the terms corresponding to further multi-particle states. The summa- 
tions run over all distinct solutions of the Bethe-Yang equations 

Qk; ai ...a n (L\0i, ...,0 n ) = m ak Lsmh8 k + ^ -i\ogS akai (9k - 9i) = 2nlk , 4eZ 

l^k 

that have total momentum zero. Using the results of [4] one can write 



3 , . 
: expz-y?(0,0) : 

: expz-</j(0,0) : 



Bi(0) 
B 3 (0) 



FlV2(0) + O (e-^) 



f» +Q(e -^) 



: expz-y>(0,0) 



B 1 (9 1 )B 2 (9 2 ] 



fH 2 (6i,6 2 ) 



+ O (e-" L ) 



8 , s 
■ expz-^(0,0) : 



Bi^B^Bi^) 



Vpi2(L\ei,e 2 ) 
Vpm(L\e u e 2 ,e 3 ) v ; 



The density factors p are obtained as 

Ph...i n {L\0\, ...,0 n ) = det 



dQ k ;ai ...a^ 

89, 



k,l=l,...,n 



In particular, for the one-particle densities we obtain 

Pk(L\0) = m^L cosh 9 



(2.2) 



The next step is to take the limit L — > 00: the exponential corrections can be dropped and 
the summations substituted with integrals 



(10 



—00 
00 



2tt 



-~Pi2(L\9i) 



£ - /If 

01,02 °° 
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where p denotes the density of zero-momentum states. For the first integral, it can be obtained 
by inspecting the Bethe-Yang equations 



miLsinh#i — Hog S^Oi — 62) 
m 2 £ sinh 2 - i log Si 2 (02 - #i ) 



2vr/i 
2irl 2 



with S12 denoting the -B1-B2 S-matrix (cf. eqn. (|A.4|) ). The second equation is actually 
superfluous due to the zero-momentum constraint mi sinh 0i + 777,2 sinh 02 = 0. Taking the 
derivative of the first equation gives 



Pi2(L\0i 



miLcosh0i + ( 1 -| — - I $12 1'/ 1 



mi cosh 0i 







^12(0) = -^ log Si 2 (0) 



using the zero-momentum constraint during the differentiation. On the other hand, the density 
factor pi2 is 

pi2(L\6i, 62) = miL cosh 0im2-C/ cosh 02 + (miLcosh0i + 1712L cosh 6*2)^12(^1 — 02) (2.3) 
and therefore 



Pu(L\t 



1 



A similar calculation yields 



pi 2 (i|0i,02) m 2 Lcosh02 
Pm(^|0i,02) 1 



pin(L|0i,02,0 3 ) miLcosh0 3 



The end result is that the correction is proportional to the volume L, and therefore it represents 
a correction to the bulk energy density 



5£ 



SE (L) 



f} /2 (0) F* /a {0) 



L 



mt 



+ 



+ 



+ 



r°° dth 
J -00 2^ 



F, 



U 2 ( 



+ 



h,02 



-,1/2, 



m 



1 r°° r°° d$! 

3! i-oo J -00 2^ 
...}+0(A 4 ) 



(mi cosh 0i + m,2 cosh 02)777,2 cosh 02 

/ 



d02 
2vr 



I'nfi 11 .- 1 ' ";; 



02=— arsinh(mi sinh 61/7712) 
2 

) 



(mi cosh 0i + mi cosh 02 + mi cosh 3 )mi cosh 03 



(2.4) 



03=— arsinh(sinh $i+sinh 9%) 



where the form factor functions are defined in Appendix [A] (the combinatorial factor in the 
last integral takes into account that states that only differ in the ordering of the rapidities are 
eventually identical). 
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R = y/4^/f3 


1.5 


1.9 


2.2 


2.5 


b 2 (TCSA) 


0.81 


1.53 


2.21 


3.04 


b 2 (FFPT) 


0.82 


1.53 


2.20 


3.01 


b sS (FFPT) 


0.0025 


4 x 10~ 5 


1 x 10~ e 


1 x 10~ 8 



Table 2.1: Comparing vacuum energy density from FFPT to TCSA numerics. The parameter 
R is related to the compactification radius of the ultraviolet limiting c = 1 free boson conformal 
field theory. 



The bulk energy density corrections can now be evaluated explicitly in units of the soliton 
mass M. Introducing also the dimensionless coupling [12] 

t = \M- 2+ P 2 / w * 2 

we can write 

The results of second order FFPT are summarized and compared to numerical values ex- 
tracted from TCSA in Table 12.11 The accuracy of the data in the table corresponds to the 
estimated precision of the TCSA results; at this level, the contribution of the integral terms is 
negligible. The deviation between FFPT and TCSA comes from two sources. For lower values 
of R, TCSA was observed to converge slower, thereby limiting the accuracy of the numerical 
determination. Albeit there exists a renormalization group method for improving convergence 
[THUS], implementing it comes with a cost (in terms of programming and running). It also 
does not seem to gain much compared to the simple-minded approach of evaluating bulk en- 
ergy by the simpler method which was applied with success in many previous examples [9lfT2]. 
Our method (also used in [H [12]) is to find the scaling regime where the ground state level 
is most linear (the region where its second derivative in L is smallest) and evaluate the slope 
of the line there. Similarly, masses can be evaluated in the region where the gap between 
the appropriate excited state and the ground state becomes closest to constant (found by 
searching for the minimum of the first derivative) and taking the value of the gap there as the 
approximate mass. 

For higher values of R, the spectrum of the theory becomes more and more dense as the 
sine-Gordon model is increasingly attractive (at the point R = 1.5 there are three breather 
states in the spectrum, while at R = 2.5 there are already eleven of them), therefore there 
are more multi-breather states to be included, and in addition there are also states containing 
solitons. 

To demonstrate that solitons contribute very little, let us also compute the value of the 
first solitonic correction, which comes from the soliton-antisoliton two-particle state. It can 
be written in a form very similar to the B\B 2 term: 

_ A 2 r de\F^ 2 (9,-e)-F s -, 1/2 (e,-9)\ 2 /4 



oo 



27r 2M cosh 9 M cosh I 



where F^ 1 ^ 2 is given in (IA.5jl . The contribution of this integral to b 2 is denoted b sS and is 
shown separately in table 12.11 The reason for the smallness of this integral is that it has a 
very limited effective support. The integrand is eventually symmetric in 9, and form factors 
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generally vanish on threshold (9 = 0). On the other hand, the form factor combination in the 
numerator exhibits an exponential decay for large 6 

\Fll 2 {6, -9) - F^ /2 (9, -6)\ 2 ~ exp 

where £ < 1 in the attractive regime. Together with the denominator this makes the integrand 
decay very fast with increasing 9. Similar behaviour happens in terms with larger number of 
particles, ensuring the convergence of all multi-particle integrals involved. Similar arguments 
hold also for the B1B2 term, but that is made larger by the appearance of smaller masses 
(mi, 1712 instead of M) in the denominator. 

The issue of whether the summation over the states with increasing number of particles 
implied by (12. lj) converges is more subtle since it is also necessary to take into account the 
various numerical prefactors (form factor normalization etc.) entering the individual contribu- 
tions and is not considered here in detail. Just as in the above calculation, the contributions 
can be naturally ordered by the sum of the masses of the constituent particles in the inter- 
mediate state, and explicit numerical evaluations support the observation that they decrease 
very rapidly when going to more and more massive states. 



3 Mass correction 



Let us now turn to evaluating the mass correction for the first breather B\. In finite volume, 
the B\ one-particle state is just the next energy level |1) above the vacuum |0)in the zero- 
momentum, zero topological charge sector. Therefore 



SE 



: exp if y?(0, 0) : 



E 



(o) 



E 



(0) 



(3.1) 



J k ^1 

and the correction to the mass gap is obtained by taking the difference to the vacuum level: 

6mi = lim 8Ei(L) - 5E (L) 

L^oo 



3.1 Bulk contributions: a puzzle and its solution 

In particular, terms linear in the volume are expected to cancel, leaving us with a finite 
correction to the mass gap. However, right with the first term a serious problem appears. The 
first contribution to (13. ip is given by the vacuum state and can be written as 



X 2 L 2 



Bi(0) : exp if 99(0,0) : 







mipi(L\0) 



+ 0(e 



X 2 L 



Fl /2 (0) 



(3.2) 



which has the wrong sign to cancel the corresponding contribution to the vacuum energy 
density i.e. the first term in eqn. (|2.4I) . 

The puzzle can be solved by observing that the contribution from B\B\ two-particle states 
(which are naively of order L° for large L) diverges as L — ► oodue to a disconnected piece. 
Such divergent pieces are finite for L < 00 , and have a dependence of L to the power of 
the number of particles involved in the disconnected part. In this case it leads to a piece 
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proportional to L, and we proceed to show that it gives the correct contributions to account 
for the mismatch noted above. The corresponding term can be written as 



9 



Bi(0) : expif<£>(0,0) : Bi(0)J3i(- 



2m\ cosh 9 — m\ 
and using the results of [4] this can be written as 



+ O (e-» L ) 





Fig {in, 9,-0) 


2 


Pi(L\0)pn 


{L\9,-9)(2m 1 cosh 9- mi) 



O (e-^) 



(3-3) 



with pi as in (|2.2p and 



pu(L\9i, 9%) = m 1 L cosh 9\ cosh 02 + mi-L(cosh0i + cosh 02) c &i 
*n(0) = -i^logSii(0) 



where 



5n(6») 



sinh + i sin 7r£ 



(3.4) 



sinh 9 — i sin tt^ 

is the -Bi-Bi scattering amplitude. A simple calculation similar to that in the previous subsec- 
tion gives the density of zero total momentum states as 



pn(L\9) = miLcoshfl + 2$n(20) 



Pii(L\9, 



m\L cosh 9 



Naive application of the infinite volume limit to (13. 3p gives 



d9 



< 2 (ivr,0, 



27T mf cosh 9(2 cosh — 1) 

However, the integral is divergent due to kinematical poles of the form factor at = 0. The 
form factors in an integrable quantum field theory satisfy a number of axioms (for the details 
we refer to Smirnov's review |13j). among which there is the kinematical residue axiom of the 
fornil 



i Res F. 



n+2( 



+ ITT, 9 ,9\,..., 9 n )iji 1 ,__i ri 



which results in the following singularity 



< 2 ("T,0, 



l-5 ij \[Su k {e-9 k ) F n c 



o. 



16 



Fi /2 (0) 



9°- 



+ 0(9°) 



(3.5) 



2 This form of the axiom is valid for self-conjugate particles; for charged particles it involves the charge 
conjugation matrix. 
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using the fact that S'ii(O) = — 1 which expresses the Pauli exclusion principle satisfied by the 
B\ particles. Therefore one must return to a more careful evaluation of the sum (13. 3p . The 
quantization of in a finite volume is given by 



m\L sinh 9 + 5 U (2(9) = 2vr J 



(3.6) 



where the quantum number is shifted by —1/2 due to the following identification of the two- 
particle phase-shift: 



S n (9) 



As a result we obtain that for fixed J 



and the leading term in the sum (|3.3p can be written as 

2 



Fl /2 {0) 



Using the identity 



we obtain 



/ miL\' 

m\L{m\V) 2 (2m\ — m\) \ 2irJ J 

v 1 - 71-2 



JGN+1/2 

2 

—X 2 L- 



fI /2 (o) 



which exactly compensates for the mismatch caused by the "wrong sign" in eqn. (|3.2|) . The 
correction to the mass term comes from the subleading L° term in the sum (|3.3p . which 
requires a very careful evaluation that is carried out in Appendix [Bj 

Moving to the next correction (-B3 term) to the bulk energy density in (12. 4ft and keeping in 
mind the above example, it is easy to see that its counterpart arises from the B\B^ contribution 
to (|3.ip . Here we encounter a different mechanism for the generation of the bulk term. The 
appropriate sum to evaluate is 

2 



61 



5i(0) :expifp(0,0) : B 1 (e 1 )B 3 (e 2 ] 



(mi cosh 6\ + 771,3 cosh 62) — mi 

with mi sinh 6\ + m.3 sinh 82 = 
It turns out that due to Si3(0) = +1 the form factor 



(3.7) 



? l/2 
113 



(iir,< 



is regular as 9\ — ► (and therefore also 62 ~ miQ\/m$ 
converts directly to an integral of the form 

dOi 



X 1 



Fl( 2 {iir,e 1 ,92) 



0) and so the above discrete sum 
\ 



2tt I (mi cosh 6*i + m3 cosh 6*2 — mi)m3 cosh ^2 



7 



(3.8) 



62=— arsinh(mi sinh di/m^) 



9 



However, the Bethe-Yang equations 



mi-Lsinh^i — i log Sis {6% — 62) = 2ttI\ 
m3-Lsinh02 - zlogSi3(#2 — 0\) = 2v:I 2 

have the solution 9\ = 62 = for I\ = I2 = 0, which is allowed due to £13(0) 
results of the paper [5], the finite volume matrix element can be written as 



fli(0) 



expz-</?(0,0) : 



fli(0)B 3 (0) 



1 



Vpi3(^|0,0)pi(L|0) 



+ 1. Using the 



F 1 1 1 / 3 2 (ivr,0,0)+m 1 LF 3 1/2 (0) 



with pi 3 obtained from (|2.3|) by replacing the index 2 with 3. The Q\ 
therefore takes the form 



term of (l3~Dl 



F 3 1/2 (0) 



^Ko,oK ,z (o) 



1/2. 



-L + 25Re 



(mi + m 3 ) 



F 3 1/2 (0) 



$13(0) 



rn-im3 



mirrii 



+ 0(L~ 







$13(0) = log Sl 3 (0) 



(3.9) 



The first term is just the correct bulk contribution, the next two terms must be added to the 
mass correction and the L~ l corrections can be discarded. 

The other bulk terms in (12. 4p . written as integrals, are expected to follow from terms of 
(|3.1|) with B\B\B2 and B\B\BiBi as intermediate state. However, their evaluation is rather 
tedious and will not be pursued here. The corresponding bulk parts in (12. 4p are small and 
therefore it is plausible that their contributions to the mass shift are small as well, compa- 
rable to numerical accuracy of TCSA and the errors made by neglecting other states. This 
assumption will be justified by the later comparison to TCSA. 



3.2 Evaluating the mass correction 

Using the formulae in the previous subsection and the end result (|B.2|) of Appendix [Bj the 
correction to the first breather mass can be written as follows (the particle composition of the 
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R = y/4^/f3 


1.6 


1.9 


2.2 


2.5 


a 2 (TCSA) 


3.8 ±0.3 


4.7 ±0.2 


6.1 ±0.1 


7.6 ±0.1 


a 2 (FFPT) 


3.66 


4.91 


6.23 


7.82 



Table 3.1: Comparing the mass correction coefficient a 2 from FFPT to TCSA numerics. The 
parameter R is related to the compactification radius of the ultraviolet limiting c = 1 free 
boson conformal field theory. The values of R are chosen to lie in a range to ensure a sufficient 
precision for the TCSA determination, for which an estimate of the numerical uncertainty is 
shown. FFPT values from (|3.10l) are reported with two decimal places accuracy. 



contributing intermediate state is indicated below each term): 
5mi 



5m\ 



(ii) 



5m\ 



(2) 



5m\ 



(i) 



5m 



(13) 



5m\ 



(22) 



5m^ 1 ^ + 5m^ + 5m^ 3 ^ + 5m^ 22 ^ + 5mf^ + 



BiBi 



B 2 

00 d9 
2k 



( 



B1-B3 



B 2 B 2 B 4 
2 



16 



A 2 x 16 



f} /2 (0) 



Fl( 2 (iTT,0) 



mf cosh 9(2 cosh 9 — 1) m\ sinh 2 cosh 9 
2 /$ u (0) 1 



4mf Am\ 



m\m 2 (m 2 — mi) 

2 



Z7T,0) 



m\miSjrit± — m\) 



m\m% 



m\m\ 



Mi 
~2^ 



( 



Fll 2 (iTr,9 1 ,8 2 ] 



\ 



(m\ cosh Q\ ± m-3 cosh 9 2 — m\)m?, cosh 9 2 



J 



(3.10) 



9 2 =— arsinh(mi sinh 61/7713) 



2 1 d6> 

2! i-oo^ 



FH 2 (iir,9,-9) 



\ 



(2m 2 cosh — mi)m 2 cosh 



/ 



(the evaluation of the terms 5m\ , (5m^ and proceeds by the already discussed meth- 

ods; there are no disconnected pieces in any of them). The mass correction can be parametrized 
with the dimensionless coefficient a 2 defined by 



5 m i 



-a 2 t 2 + 0(t 4 



which is compared to TCSA data in Table [371 
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4 Conclusions 



In this paper it was shown how to use finite volume techniques to go beyond first order in form 
factor perturbation theory. The second-order corrections to vacuum energy and first breather 
mass was evaluated in double sine-Gordon theory. In principle, the method works for higher 
corrections, and for other quantities, such as the S matrbH as well. 

The results of second order FFPT are in good agreement with numerical data from TCSA. 
In addition, the regularization techniques developed to evaluate disconnected contributions 
can also be used for the form factor expansion of finite temperature two-point correlators. 
Eventually, during the typing of this manuscript there appeared an independent work by 
Essler and Konik [16], which uses similar finite volume techniques for correlators, and also 
introduces another, novel infinite volume regularization procedure. 
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A Breather form factors in sine-Gordon theory 



To obtain matrix elements containing the first breather, one can analytically continue the 
form factors of sinh-Gordon theory obtained in [17] to imaginary values of the couplings. For 
the theory obtained by setting A = in (|1.2jl . the result reads 



'!)•••) Vn 



where £ = /3 2 /(8vr - /3 2 ), 

Qa ^ { x l ) • • • > x n) 

HO 



g a {p) [a)s(f\(0) n U 



B 1 (fi 1 )...B 1 (0, 

ft(0j-(h) 



i<j 



. e Oi _|_ QOj 



det [a + i- j}{ ofi-j^ ■ ■ ■ . x ™) ij= i,..., n _i if n > 2 
sin 7r£a 



sin 7r£ 



< „ . < 

2 cos — A / 2 sm — exp 



dt t 



2tt sin t 



(A.l) 



and 



/ e (0) = v{m + 9, -l)v(in + 0, -£)v(in + 9,1 + £)v{-m - 9, -l)v{-m - 9, -^)v(-iir - 9, 1 



<ar\ fr ( e + m{2k + Q \' 



k=l 

x exp ( r * ( < - — + (N + 1 - Ne~ 2t ) e~ 2 ^ 

\Jo t V 4sinh| 2vrcoshf v ; 2sinh 2 i 



3 The evaluation of S matrix corrections can be carried out by calculating the shifts of two-particle levels, 
which depend on the phase shift in finite volume. 
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gives the minimal B\B\ form factoiQ, while denotes the elementary symmetric polynomial 
of n variables and order k defined by 



Furthermore 



i=l 



k=0 



9aiP) = (e ia ^) 



x exp 




2sinh ( ^t) cosh ( ( 1 - ~ ) / ) sinh / 



sinh 2 (£t) 



4tt 



is the exact vacuum expectation value of the exponential field [18], with M denoting the 
soliton mass related to the coupling \i via [T9] 



2r(A) 



i 

2-2A 



M 



2-2A 



13 2 
8vr 



(A.2) 



Formula IjA.lj) also coincides with the result given in |20J. 

Form factors containing higher breathers can be obtained using that B n is a bound state 
of B\ and B n _i; therefore sequentially fusing n adjacent first breathers gives B n . Following 
the lines of reasoning of Appendix A of the paper [21] one obtains 



Ft 



ki...krnh...l s \ u l-> ■ ■ ■ } a r, v, 



(A.3) 



Jaf3<S>(0) 



B kl (0i) ■ ■ ■ B kr (9 r )B n (9)B h (9[) . . . B ls (9' s ) 



7?i7l 3 2 • • • lln-l 



fcl...fcrll...lii...i s I 01 1 • • • > "r,v + 



1 — n 



3 — n 



iir£,...,9 + 



n — 1 



) * * * ? s 



where 



is the B\B k 



I lk 



2 tan ^p- tan 



(fc+l)7r£ 



tan 



B k+ i coupling, defined as the residue of the appropriate scattering amplitude: 

2 



Res S lfc (0) 

fl_ i7T(fc + l)g 

(7— 2 

sinh 9 + i sin Ek^f!ii» Q sm 1 ^ 

sinh — i sin Ek^ijj^ smh 9 — i sin 1 ^ 



(A.4) 



4 The formula for the function v is in fact independent of N; choosing N large extends the width of the 
strip where the integral converges and also speeds up convergence. 
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Using the results of [20] we also quote here the simplest solitonic form factor, which is needed 
in the main text 



, e+j-rr 



(A.5) 



where 



G(0) = iC\ sinh - exp 



At sinh^(l - f )sinhi(£ - 1) 
t sinh2i cosht sinhi£ 



Cl 



exp 



dt sinh 2 \ sinh t(£ - 1) 



Iq t sinh2i cosht sinht^ / 

B Evaluation of disconnected contributions 

To obtain #774 we need to evaluate the 0(L°) of the sum (13,3(1 . which is 







2 


Pi(L\0)pu 


{L\9,-e)(2m 1 cosh 9- mi) 



+ O (e-" L ) 



Using 



16 



F x 1/2 (0) 



e 2 



+ 0(6°) 



we can subtract the singular piece to write 
/ 



-A 2 L 2 £ 





Flll{m,e-6) 


2 


Pi(i|0)pn 


{L\0,-e)(2mico$b.O - mx) 



16 



F! 1/2 (0) 



which can be readily converted for L — > oo to the following integral: 



f 



d6 



16 



Fl /2 {0) 



2ir I mf cosh 6* (2 cosh 6* — 1) m\ sinh 2 cosh 9 
Therefore what we need is the 0(L°) part of the subtracted term 



16 



Fi /2 (0) 



/sinh J 0pi(L|O)pii(L|0,-0)mi 



(B.l) 



where 



Pi(L\0) = m x L, p u (L\e,-9) =m 2 L 2 cosh 2 (9 + 2miLcosh(9$ii(2(9) 
According to eqn. (|3.6|) the rapidity is quantized as 

miLsinhfl + £n(20) = 2vrJ 
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where J is a positive half-integer. Using the arguments of Subsection 13.11 



16 I Fi /2 (0) 



16|i?i 1/2 (0) 



y2j2 " I 1 ( m lL \ _ ^2j2 " 1 1 v ' y 

^ mfL(miL) 2 \ 2-7rJ / ^ m\L(miL) 



1 



sinh + 



X 2 L 



fI /2 (o) 



is just the O(L) part of the sum (|B.lj) . which can be explicitly subtracted without affecting 
the O(L ) part. Thus the O(L ) term of (jRljl can be obtained as the L — ► oo limit of 



16- 



^i 1/2 (o)| 



16- 



F x 1/2 (0)| 



E 
E 



sinh 2 cosh (m-iL cosh 9 + 2$ii(26>)) miLsinh 2 



1 



V 



sinh 9 



where an intermediate subtraction was inserted to simplify the evaluation. The second sum 
can be written 



E 



i 



i 



V 



sinh 2 9 



sinh + M2*)V 

rn\L J 



Em 2 L 2 



26n(20) 



Air 2 J 2 \ miLsinh( 



+ 



Mgg) 

mi L sinh 9 



Now the singular part is in the J -2 prefactor, and the remainder is finite for any fixed J when 
L — ► oo. Since in the infinite volume limit # — > for any fixed J, and the summation over J 
is uniformly convergent, permitting to exchange the limit with the sum, we can write 



E 



i 



i 



sinh 6 

V 

For the first sum we obtain 

E 



E 5j5 4 *»<°> + °( i0 » = ^^^Y^l + 0(L") 



1 



sinh 2 9 cosh 9 [m x L cosh 9 + 2$ u (20)) mi L sinh 9 



Y — 

y> sinh 2 9 
1 $n(0) 



4 4 
where in the first term we used 



1 - 1/ cosh 2 9 
+ OiL- 1 ) 



+ 2- 



gii(2g) 
2 L 2 cosh 3 



+ 0(L 



-3 



E 



1 - 1/ cosh 2 (9 
mi-L sinh 2 



o 



oo jn 

— (m 1 Lcosh^ + $n(26»)) 

27T 



= I + o^- 1 ) 

mi L cosh 4 



(the integrand is non-singular at 9 = 0, hence it can be evaluated by a density integral), while 
for the second term: 



E— 



2*ii (20) 



sinh 2 9 m 2 L 2 cosh 3 



E 



2$n(0) 
m 2 L 2 sinh 2 9 
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where again, all parts non-singular at 6 = were evaluated at the origin. Collecting all the 
pieces we obtain that the subtracted part (jB.ip equals 



\ 2 L 2 Y- 

Z 4 



16 



Fl /2 (0) 



sinh 2 9 p 1 (L\0)p 11 {L\9,-9)m 1 
and therefore 



X 2 L ■ 



Fi /2 (0) 



-+A i xl6 



rat 



Fi /2 (0) 



$n(0) 1 
4m^ 4mf 



5m 



(ii) 



-A : 



Jo 



de 

2^ 



( 



FlS(in,9,-9) 



16 



-,1/2 



(0) 



-A 2 x 16 



V 

F x 1/2 (0) 



mf cosh 0(2 cosh — 1) mf sinh 2 cosh 6* 
2 /$ n (0) 1 



(B.2) 



4mf 



4mf 
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